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Abstract 

The spatially-temporal model of an all-solid-state passively Q-switched os- 
cillator with an active medium providing the stimulated Raman scattering is 
presented. The model does not presume a Gaussian shape of the cylindrically 
symmetric modes at both fundamental and Stokes wavelengths. It is found, 
that the highly-nontrivial spatially-temporal dynamics can be regularized by 
the optimal choice of the oscillator parameters, viz. initial transmission of a 
saturable absorber, curvature of a spherical mirror, and output mirror trans- 
mission at the fundamental and Stokes wavelengths. As a result, the pulse 
can be substantially temporally squeezed and spatially broadened at both 
fundamental and Stokes wavelengths. 
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1. Introduction 

Solid-state Q-switched oscillators allowing nano- and sub-nanosecond puls- 
ing find applications in a lot of areas including medicine, spectroscopy, en- 
vironment monitoring, etc. Passive Q-switching based on the use of both 
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semiconductor [l|, ^ and crystalline saturable absorbers is of particular 
interest due to compactness, simplicity, high damage threshold, and diode- 
pumping ability of an oscillator. 

Among the active media allowing diode-pumped Q-switching, KY(W04)2 
(KYW) and KGd(W04)2 (KGW) crystals doped by Yh^+ and ions are 
known as the materials providing an efficient stimulated Raman scattering 
(SRS) 0, i H H 8- As a result, there is possible a high-efficient self- 
frequency shift and a simultaneous two-wavelengths pulsing (e.g., at 1.35 
and 1.54 /im for a Nd:KGW active medium) directly from an oscillator. 

Theoretical studies of a Q-switched oscillator with the intra-cavity SRS 
have been based on the well-established rate-equations approach (ol. [lol. 11 . 
The oscillator parameters providing a pulse width minimization [lO(| and 



an output energy maximization [llj at both fundamental and Stokes wave- 
lengths have been defined. Simultaneously, it has been found that the spa- 
tial structure of a laser field is strongly affected by the SRS jof. This means 
that the mode transformation has to be taken into account along with the 
temporal evolution of fields and populations inside an active medium and a 
saturable absorber. Such a model has been developed to a Gaussian mode 
approximation for both fundamental and Stokes fields 12| • This has allowed 



defining the optimal values of the Raman gain and the ratio of pump and 
laser beams, which provide a single pulse operation of oscillator with the 
most efficient SRS. 

Nevertheless, a transient character of Q-switching can prevent from the 
CW mode formation and disturb substantially the spatial structure of a laser 
field. This requires to take into account the spatial dynamics on a par with 
the temporal one. 

In this work we present the analysis of spatially-temporal dynamics of 
a passively Q-switched Raman-active oscillator. The analysis is based on a 
rate-equations approach but without imposing a limitation on the shape of 
cylindrically symmetric transverse distribution of a field. The results demon- 
strate that there exist both spatial extra-broadening and squeezing as well as 
transition between non-Gaussian and Gaussian spatial profiles in dependence 
on the oscillator parameters (viz. output mirror transmission at the funda- 
mental and Stokes wavelengths, initial transmission of a saturable absorber 
and curvature of a spherical mirror). The SRS contribution controlled by 
the optimized sets of oscillator parameters allows the regularization of the 
spatially-temporal structure of a field and the substantial pulse shortening. 
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2. Numerical procedure 

The oscillator under consideration consists of the flat (output) and spher 
ical mirrors. The active medium (Nd'^+:KGW) is placed at the resonator cen- 
ter, and the saturable absorber (V^+:YAG) is placed on the output mirror. 

The model is based on a direct generalization of that presented in 



It is supposed that the time-dependent (t is the time, which is periodical with 
the cavity period Tcav- t G [0, Tcav]) fields af^s (indexes / and s correspond to 
the fundamental and Stokes fields, respectively) have the radially-symmetric 
transverse spatial distributions (r e ]0,R] is the radial coordinate). 

Inside an active medium, the dynamics is modeled on basis of the split- 
step Hankel's method. That is the active medium volume is considered to be 
divided into 10 transverse slices with the thickness Az {z is the longitudinal 
coordinate), in which the dynamics within the time-domain and the domain 
of spatial frequencies {u = c/27r_R, c is the velocity of light, R = 0.5 ^ 1.5 
cm) [l3| is evolved step- by-step: 



CLf {z, r, t) = 7^^(^f t) [g (r, t) - 5-, |a, [z, r, t) |^] 

cav 

as {z, r, t) = 7^^(^s (z, r,t)gs\af {z, r, t) |^ 

cav 

g{r,t) = -^^g{r,t) \af iz,r,t)f 



> ® (1) 



)a/,s {z + Az, z/, t) = ttf^s {z, u, t) exp ^—ikj^sAz + -kf^sAzX'j gu'^^ 



We shall suppose, that the field intensities \af^sf are normalized to he/ XgCgTcav 
[h is the Planck constant; \f =1.35 /xm and =1.54 /xm are the fundamen- 
tal and Stokes wavelengths, respectively; ag =0.76x10"^^ cm^ is the gain 
cross-section, T^av =0.8 ns), and the time is normalized to Tcav The gain 
coefficient is g and its initial value equals to the threshold one: (IuI/Tq + 
Inl/p/ + l)/2Lg. In the last expression, the varied values Tq and pf corre- 
spond, respectively, to the initial transmission of saturable absorber and the 
reflection of output mirror at the fundamental wavelength; / = 0.05 is the 
unsaturable net-loss coefficient. Lg =5 cm is the active medium length. The 
stimulated Raman scattering inside an active medium is described by the 
coefficient gg and its dimensional value amounts to 6 cm/GW. It is assumed 
that there exists no anti-Stokes and higher-order Stokes scattering as well as 
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that the gain is saturable by only fundamental field. Also, the initial gain 
distribution is assumed to be spatially homogeneous. 

The field representations aj^s within the spatial frequenc y d omain are 
obtained by the means of the fast Hankel's transformation 13|, [ij]. kf^s 
are the wave-numbers corresponding to the fundamental and Stokes fields, 
respectively. 

A similar procedure describes the dynamics inside a saturable absorber: 



d/ {z,r,t) 
hs iz,r,t) 
n (r, t) = - 



Az 



2T 



cav 

Az 



■ttf {z, r, t) [n (r, t) - aesa {N - n (r, t))] 
as {z,r,t)n{r,t) 



2T 

■^-^ cav 

■^n (r, t)\af{z,r,t)\' - -j-^n (r, t) |a, {z, r,t)f + ^^^^ 



> ® (2) 



®a/,s {z + Az, z/, t) = a/,s {z, v,t) exp y-ikf^sAz + -kf^sAzXj,,iy" 

Here, n is the loss coefficient and its initial value amounts to = (In 1/Tq )/2La 
(La=0.1 cm is the saturable absorber thickness). The loss is saturable by 
both fundamental and Stokes fields (a/ = 95, as = 4, when the normaliza- 
tion to (jg is presumed). Excited-state absorption is taken into account for 
the fundamental field [a^sa =0.1). = 26 is the loss relaxation time (the 
normalization to Tcav is presumed). 

Free propagations between the active medium and the absorber as well 
as between the active medium and the spherical mirror are considered in the 
frequency domain: 



df^s {z + z/, t) = df^s {z, t) exp \^-ikf^sL + -fc/,^LAj ,,z/^J , (3) 

where L = 5 cm is the propagation length. 

The refiection from a spherical mirror obeys: 

af^s{z,r,t) = ttf^s {z,r,t) exp [ikf^sr'^ / Rm) (4) 

where Rm is the variable radius of the mirror curvature. 

The loss on an output mirror as well as the net unsaturable loss are taken 
into account by means of the following mapping: 
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Figure 1: a: Contour plot of the output intensity at the fundamental wavelength from a 
Q-switched oscillator with To =0.3, pf —0.7, ps =0, Rm =300 cm. h: Output intensity 
integrated over r (integral power; upper plot) or t (integral fluency; lower plot). 



af^s (z, r, t) = aj^s (z, r, t) exp [-0.5(/ + ln(l/p/,,))] , (5) 

where ps is the output mirror reflection at the Stokes wavelength. 

In the simulations, the time window is approximated by 100 points and 
the window of transverse coordinate is approximated by 2000-i-6000 points 
in dependence on the R value, which varies from 0.5 to 1.5 cm in order to 
exclude the boundary effects. The solution in the time domain is evaluated 
using a fourth-order Runge-Kutta method. The initial dimensionless fields 
have the Gaussian spatial profile with the size, which equals to that of the 
fundamental CW mode at A/. The initial dimensionless amplitudes o/ and 
ttg amount to 10~^ and 10~^°, respectively. 

3. Spatially-temporal structure of laser field from a Q-switched 
oscillator 

The simulations based on the model described in the previous Section 
demonstrate that the spatially-temporal structure of laser field from a Q- 
switched oscillator depends non-trivially on the oscillator parameters, which 
effect on the field dynamics at both fundamental and Stokes wavelengths. 

Fig. [H a shows the contour plot of the output intensity at the fundamen- 
tal wavelength in absence of an effective SRS in an oscillator {ps =0). The 
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Figure 2: Contour plots of the output intensity at the fundamental (a) and Stokes (6) 
wavelengths. Tq =0.3, pf =0.7, ps =0.7, Rm =300 cm. 

intensities are averaged over the cavity period in Figs. [H a; [2], HI [5], [3, and [HI 
Both low initial transmission of a saturable absorber and comparatively small 
radius of curvature of a spherical mirror initiate the multiple pulse dynamics, 
which is clearly visible in Fig. [H h (upper plot) showing the output intensity 
integrated over the transverse spatial coordinate. The greatest pulse is com- 
paratively long (f»5 ns) and the narrow "mode" (the radius equals to ^130 
/im) has an elongate wing (lower plot). As will be shown, a non-Gaussian 
transverse shape is typical for the regime under consideration (see (ol). 

When a field at the Stokes wavelength is locked inside an oscillator owing 
to non-zero p^, the dynamics changes. For the case presented in Fig. [H the 
multipulsing disappears (Figs. [2], a and [3]) and the pulse at the fundamental 
wavelength shortens substantially (;^1.2 ns, see Fig. [3l upper plot) fioj. The 
mode at the fundamental wavelength widens and becomes Gaussian (Fig. 
El lower plot; see also j^). Simultaneously, the short pulse at the Stokes 
wavelength develops (Figs. [21 b and [3]). Such a pulse is located on the tale 
of the fundamental pulse (Fig. [3l gray curve in upper plot) and has a broad 
trapezoidal transverse profile (the mode radius reaches k,1A mm in Fig. [3l 
gray curve in lower plot). 

Thus, the transverse distribution excesses substantially the CW mode 
size, which equals to ^265 /xm for the spherical mirror under consideration. 
Such a spatial extra-broadening can be attributed to the saturable absorber 
saturation, which results in an excitation of high-order spatial frequencies. 
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r, cm 

Figure 3: Output integral power (upper plot) and integral fluency (lower plot) at the 
fundamental (black curves) and Stokes (gray curves) wavelengths. Parameters correspond 
to Fig. H 

A spatial extra-broadening is possible at the fundamental wavelength, as 
well. Fig. m a demonstrates such an extra-broadening (up to ~0.9 cm). 
The field covers almost a whole active crystal radius. It should be noted, 
that the output CW mode radius equals to ?«440 /im for the spherical mirror 
under consideration. The Q-switching "mode" has a shape of a dilative ring, 
which rises and then disappears (after ns) with a gain depletion. In the 
presence of SRS, such a "mode" is suppressed owing to an inhomogeneous 
"collapsing" spatially-temporal behavior at both fundamental and Stokes 
wavelengths (Fig. |U h and [15]). 

An effect of the SRS on the Q-switching dynamics allows controlling the 
temporal and spatial profiles of a pulse at both fundamental and Stokes 
wavelengths. An imbalance between the loss and gain saturation at the 
fundamental wavelength causing the multiple pulsing can be compensated by 
the nonhnear loss due to SRS. As a result, the well-shaped and substantially 
shortened [lo| pulses with the broad spatial profiles at both wavelengths 
appear (Figs. 

When the oscillator parameters are not optimized, the spatially-temporal 
structure of the Stokes field can be complicated. For instance, the radial 
position of the SRS spike can depend on the time (Fig. [3, a, where the 
SRS spikes move off the optical axis with the pulse evolution; this effect 
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Figure 4: Contour plots of the output intensity at the fundamental wavelength. Tq =0.3, 
Pf =0.5; =0 (a), 0.7 (6); Rm =2000 cm. 
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Figure 5: Contour plots of the output intensity at the fundamental (a, b) and Stokes (c) 
wavelengths. Tq =0.5, pf =0.9; p^ =0 (a), 0.9 (6, c); Em =1000 cm. 
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Figure 6: Output integral power (upper plot) and integral fluency (lower plot) at the fun- 
damental (black solid and dotted curves) and Stokes (gray curves) wavelengthes. ps =0 
(black solid curves) and 0.9 (black dotted and gray curves). Another parameters corre- 
spond to Fig. [5l Dotted and gray curves are vertically rescaled for convenience. 



has been observed experimentally in [16|). A similar situation is possible 
also at the fundamental wavelength (Fig. HJ a). Even through the SRS is 
well-synchronized (i.e. the SRS appears synchronously at the different radial 
positions), the mode profile is inhomogeneous, as a rule (Fig. [3, h). 

The oscillator parameters optimization allows regularizing the spatially- 
temporal dynamics at both fundamental and Stokes wavelengths. In par- 
ticular, the spatial profile of the Stokes field becomes homogeneous though 
non-Gaussian (Fig. [8]). The numerical analysis reveals, that there exist some 
optimal sets of parameters (see Table [1]), in particular the optimal curva- 
ture of a spherical mirror (:^1000 cm in the case under consideration), which 
provide the regular pulses within a broadest range of parameters (viz. pj, 
Ps, and To). The initial transmission Tq >0.5 results in a sufficiently strong 
fundamental field, that causes an efficient SRS. The reflectivity at the Stokes 
wavelengths has to be sufficiently high (>0.7) to provide an efficient conver- 
sion from the fundamental field to the Stokes one. 

The resulting pulses at both wavelength are shortened in comparison with 
a single- wavelength regime and have the smooth and broad transverse spatial 
profiles. The last property allows using a larger volume of an active crystal. 
The "mode" shape is close to the Gaussian one for the fundamental field and 
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Figure 7: Stokes intensity profiles for (a): i?M=2000 cm, pf— 0.7, ps— 0.7, and (b): 
i?Af = 1000 cm, p/==0.9, ps=0.5. To=0.5. 




Figure 8: Stokes intensity profile for Rm=WOO cm, p/=0.9, ps=0.9, and To=0.5. 
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Table 1: Oscillator parameters providing the regular spatial-temporal pulse profiles at 
both fundamental and Stokes wavelengths. 



Rm, cm 


To 


Pf 


Ps 


300 


0.3 - 0.5 


0.7- 1 


0.7-0.9 


1000 


0.3 - 0.5 


0.5 - 1 


0.7-0.9 


2000 


0.5 


0.9-1 


0.9 



to the trapezoidal one for the Stokes field. 
4. Conclusion 

The model of a Q-switched Nd:KGW oscillator with a V:YAG saturable 
absorber, which takes into account a transverse spatial field distribution at 
both fundamental and Stokes wavelengths, has been developed. The nu- 
merical analysis has demonstrated that an imbalance between the loss and 
gain saturation causes not only a multiple pulsing but an aberration of an 
oscillator "mode" . When the initial transmission of a saturable absorber is 
sufficiently low, the spatial extra-broadening of a field develops so that a 
"mode" covers almost whole volume of an active crystal. 

Manipulations with the beam size on a saturable absorber (by means of 
the change of a spherical mirror curvature), the initial transmission of an 
absorber, and the output mirror reflectivity at both fundamental and Stokes 
wavelengths allow controlling the spatially-temporal profiles of the output 
pulses. To obtain a single pulsing at two-wavelengths with a broad and 
homogeneous spatial distribution, the initial transmission of the saturable 
absorber has to be about of 0.5, the reflectivity of the output mirror at the 
fundamental wavelength has to approach 1 (it can range within 0.5 - 1 for 
the optimal curvature of the spherical mirror, which is of ~1 m in the case 
under consideration) , and the reflectivity at the Stokes wavelength has to be 
about of 0.9. The SRS causes a temporal squeezing of the pulses at both 
wavelengths, whereas the spatial proflles are greatly stretched in comparison 
with the CW-mode size. When the oscillator parameters are optimized, the 
transverse spatial profiles are near-Gaussian at the fundamental wavelength 
and trapezoidal at the Stokes wavelength. 

The obtained results are of interest for development of the passively Q- 
switched Raman-active oscillators producing the pulses with durations of 
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about of few nanoseconds and possessing the smooth transverse spatial pro- 
files covering a considerable part of an active crystal. 
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